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ABSTRACT 



Context. The triple stellar system S Vel system presents a significant infrared excess, whose origin is still being debated. A large 
infrared bow shock has been discovered using Spitzer/MIFS observations. Although it appears as a significant contributor to the 
measured IR excess, the possibility exists that a circumstellar IR excess is present around the stars of the system. 
Aims. The objective of the present VISIR and NACO observations is to identify whether one of the stars of the (5 Vel system presents 
a circumstellar photometric excess in the thermal IR domain and to quantify it. 

Methods. We observed (5 Vel using the imaging modes of the ESOA^LT instruments VISIR (in BURST mode) and NACO to resolve 
the A-B system (0.6'' separation) and obtain the photometry of each star. We also obtained one NACO photometry epoch precisely at 
the primary (annular) eclipse of 6 Vel Aa by Ab. 

Results. Our photometric measurements with NACO (2.17//m), complemented by the existing visible photometry allowed us to 
reconstruct the spectral energy distribution of the three stars. We then compared the VISIR photometry (8.6-12.8 yt/m) to the expected 
photospheric emission from the three stars at the corresponding wavelengths. 

Conclusions. We can exclude the presence of a circumstellar thermal infrared excess around 6 Vel A or B down to a few percent level. 
This supports the conclusions of Caspar et al. (2008) that the IR excess of 6 Vel has an interstellar origin, although a cold circumstellar 
disk could still be present. In addition, we derive the spectral types of the three stars Aa, Ab, and B (respectively A2IV, A4V and 
F8V), and we estimate the age of the system around 400-500 Myr. 

Key words. Stars: individual: (HD 74956, S Vel); Stars: binaries: eclipsing; Methods: observational; Techniques: high angular reso- 
lution 
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1. Introduction 

The southern star S Vel (HD 74956, HIP 41913, GJ 321.3, GJ 
9278) is the 43^^ brightest star in the visible sky, with my = 1.96, 
and is located at a distance of 24 Ape (n = 40.90 ± 0.38 mas; 
ESA I1997I) . Amazingly, S Vel was discovered only very recently 
by Otero et al. (2000) to host the brightest eclipsing binary sys- 
tem in the southern sky, one of the very few observable with 
the unaided eye. Its orbital period of 45 days is exceptionally 
long in terms of probability of occurrence, especially for such 
a nearby star. A historical note on the discovery and recent de- 
velopments in the study of S Vel can be found in Argyle, Alzner 
& Horch (12002b , together with the orbital parameters of the A- 
B system (orbital period of ^ 142 years). In the following, we 
refer to the eclipsing pair as 6 Vel A, and to the fainter, visual 
component as ^ Vel B. The two eclipsing components of A are 
referred to as Aa and Ab, by order of decreasing brightness. 
Kellerer et al. (120071 ) excluded the physical association of the 
angularly nearby pair of stars sometimes labeled S Vel C and D 
with the S Vel A-B system. 
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An infrared (IR) excess associated with the ^Vel A-B sys- 
tem had been detected by IRAS (1986) . This indicated that S Vel 
could belong to the "Vega-like" class of objects, i.e. a main se- 
quence star surrounded by a optically thin debris disk. However, 
Caspar et al. ( 2008 ) recently discovered a spectacular infrared 
(IR) bow shock around S Vel at 24 and 70 jim using Spitzer/MIFS 
images. This very large structure (^ V) is explained by these 
authors as the result of the heating and compression of the in- 
terstellar medium by the photons from S Vel, as the trio moves 
through the interstellar medium (ISM). From a detailed model- 
ing of the star-ISM interaction, Gaspar et al. conclude that the 
bow shock contribution is sufficient to explain the observed IR 
excess of S Vel without resorting to a circumstellar debris disk. 
They also note that accretion of interstellar material (A Bootis 
phenomenon) could nevertheless take place on S Vel, although 
at a slow rate. 



In this article, we present the results of our search for circum- 
stellar IR excess in the inner S Vel A-B system. Sect. [2] is ded- 
icated to the description of our high angular resolution VISIR 
and NACO observations, whose analysis is presented in Sect. [3] 
In Sect, m we compute the magnitudes of each of the three stars 
of the system, that are used to derive their physical parameters. 
We present in Sect.[5]our estimates of the thermal IR excesses of 
6 Vel Aab and B and a short analysis of the evolutionary state of 
the system. 
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2. Observations and data reduction 

2.1. VISIR 

For the thermal-IR part of our program, we used the VISIR in- 
strument (Lagage et al. 12004b . installed at the Cassegrain focus 
of the Melipal telescope (UTS) of the ESO/Very Large Telescope 
(Paranal, Chile). VISIR is a mid-IR imager, that also provides a 
slit spectrometer. VISIR can in principle reach a very high spa- 
tial resolving power, thanks to the 8 m diameter of the telescope. 
However, under standard conditions at Paranal (median seeing 
of 0.8'' at 0.5//m), the 8 m telescope is not diffraction limited 
in the MIR (seeing ^ 0.4'' vs. 0.3" diffraction). A few mov- 
ing speckles and tip-tilt usually degrade the quality of the image 
(see e.g. Tokovinin, Sarazin & Smette 2007 ). To overcome this 
limitation, a specific mode of the instrument, the BURST mode 
(Doucet et al."2007al.l2007bl), give the possibility to record sep- 
arately a large number (several tens of thousand) very short ex- 
posures (A^ < 50 ms), in order to freeze the turbulence. The data 
processing procedure we developed to reduce the resulting data 
cubes is described in Kervella & Domiciano de Souza ( 20071) . 
During the processing, the frames were precisely co-aligned on 
the position of S Vel A. 

We observed ^ Vel and its three calibrators, HD 82668, HD 
67582 and HD 80007 in visitor mode at Paranal during the night 
of 23-24 April 2008. At that time ^ Vel was out of ecplipse. The 
sequence of the observations is presented in Table[T] We adopted 
standard chopping and nodding offsets of 8", with respective 
periods of 4 and 90s. The calibrators were selected from the 
Cohen et al. (1999) catalog of spectrophotometric standards for 
IR wavelengths, except HD 80007 (for the Nell filter). For each 
observation of S Vel and its calibrators, we selected during the 
data processing ^ 40% of the total number of frames, reject- 
ing those with the lowest Strehl ratio (estimated from the peak 
intensity in the frame). For instance, the first two observations 
#A and B of S Vel (20 000 selected frames) correspond to a total 
of 50 000 frames before selection. The number of frames listed 
in Table [T] corresponds to the result of this selection. The cubes 
were then averaged to obtain diffraction-limited images in our 
three filters: PAHl, PAH2, and Nell, that have respective cen- 
tral wavelengths of A = 8.59, 11.25 and 12.81yum. The resulting 
images of ^ Vel and the calibrators are presented in Fig. [T] The 
observation #C of HD 82668 was affected by saturation of the 
array, and was not used for our anaysis. 

2.2. NACO 

We observed S Vel at ten epochs in April-May 2008 using the 
Nasmyth Adaptive Optics System (NAOS, Rousset et al. 2003 ) 
of the Very Large Telescope (VLT), coupled to the CONIC A in- 
frared camera (Lenzen et al. ll998l) , abbreviated as NACO. These 
observations were obtained to provide high-precision differential 
astrometry of the eclipsing system ^ Vel A relative to B, and we 
will analyze them along this line in a forthcoming article. 

We selected the smallest available pixel scale of 13.26 ± 
0.03 mas/pix (Masciadri et al. 2003 ), giving a field of view of 
13.6"xl3.6". Due to the brightness of ^Vel, we employed a 
narrow-band filter at a wavelength of 2.166 ± 0.023 yum (here- 
after abbreviated as 2.17) together with a neutral density filter 
(labeled "ND2 .short"), with a transmission of about 1.5%. 

Table [2] gives the list of the observations. Each of our ten 
epochs consisted in approximately 50 short exposures grouped 
over less than 10 minutes, with a detector integration time (DIT) 
of 0.8 s each to avoid saturation. Four such short exposures were 



Table 1. Log of the observations of S Vel and its calibrators, HD 
82668, HD 67582 and HD 80007 with VISIR. MJD is the mod- 
ified Julian date of the middle of the exposures on the target, 
minus 54 580. The Detector Integration Time (DIT) is given in 
milliseconds for one BURST image. 6 is the seeing in the visi- 
ble (A = 0.5 Jim) as measured by the observatory DIMM sensor, 
in arcseconds. N exp. is the number of selected and averaged 
exposures. The airmass (AM) is the average airmass of the ob- 
servation. 



# 


MJD* 


Star 


Filter 


DIT 


N exp. 


on 


AM 


A 


0.0134 


(5 Vel 


PAHl 


20 


10000 


0.86 


1.18 


B 


0.0231 


(5 Vel 


PAHl 


20 


10000 


1.00 


1.20 


C 


0.0374 


HD 82668 


PAHl 


20 


5 000 


1.10 


1.20 


D 


0.0603 


HD 67582 


PAHl 


16 


10000 


0.84 


1.30 


E 


0.0774 


(5 Vel 


PAHl 


16 


8 000 


0.70 


1.33 


F 


0.0874 


(5 Vel 


PAHl 


16 


8 000 


0.67 


1.37 


G 


0.0977 


(5 Vel 


PAH2 


8 


10000 


0.60 


1.42 


H 


0.1093 


(5 Vel 


PAH2 


8 


12 000 


0.53 


1.48 


I 


0.1206 


(5 Vel 


PAH2 


8 


10000 


0.80 


1.54 


J 


0.1363 


HD 67582 


PAH2 


8 


10000 


1.15 


1.92 


K 


0.1480 


HD 67582 


PAH2 


8 


10000 


1.50 


2.09 


L 


0.1632 


(5 Vel 


Nell 


16 


10000 


1.29 


1.91 


M 


0.1860 


HD 80007 


Nell 


16 


10000 


1.12 


1.46 


N 


0.1979 


HD 80007 


Nell 


16 


10000 


1.00 


1.52 
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Fig. 1. Top row: VISIR average images of S Vel in the PAHl (#A, 
B, E and F in Table©, PAH2 (#G, H, I) and Nell (#L) bands. 
Bottom row: Calibrator images in the PAHl (HD 67582), PAH2 
(HD 67582) and Nell (HD 80007) bands. 



co-averaged providing 3.2 s integrated exposure time per image. 
A few images (9 over 510 in total) were affected by the opening 
of the adaptive optics loop, and were removed from the process- 
ing. We obtained nine epochs outside of the eclipses, and one 
precisely at the phase of the primary eclipse (on 18 May 2008), 
only 23 minutes before the center of the eclipse. As a remark, the 
primary eclipse (the deeper in photometry) is when the smaller, 
cooler star (Ab) passes in front of the bigger, hotter star (Aa), 
causing an annular eclipse. The secondary eclipse (shallower in 
photometry) is when the cooler star is totally eclipsed by the 
hotter one. The epochs of the primary and secondary eclipses of 
S Vel A are given by the ephemeris of OtercQ: 

Mini = HJD 2452798.557 -h 45.1501 E (1) 

Minll = HJD 2452818.200 -h 45.1501 E (2) 
^ http ://ar. geocities . com/varsao/delta_ Velorum.htm 
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Table 2. Log of the observations of ^Vel with NACO in the 
NB2.17 filter. MJD* is the average modified Julian date minus 
54000, HJD* the heliocentric JuHan date minus 2454000 and 
is the phase of the eclipsing binary (see text). N is the num- 
ber of frames. As in Table [U 6 is the seeing in the visible and 
AM is the airmass. The listed figures are median values over the 
observations. 



# 


MJD* 


HJD* 


0(MinI) 


N exp. 


on 


AM 


1 


557.0224 


557.5245 


0.9582 


51 


0.80 


1 


157 


2 


560.9976 


561.4996 


0.0463 


46 


0.91 


1 


164 


3 


562.0121 


562.5141 


0.0687 


51 


0.67 


1 


156 


4 


563.0048 


563.5067 


0.0907 


51 


0.89 


1 


158 


5 


576.9715 


577.4732 


0.4000 


51 


0.66 


1 


156 


6 


579.0231 


579.5247 


0.4455 


50 


1.06 


1 


191 


7 


580.9917 


581.4932 


0.4891 


51 


1.24 


1 


164 


8 


591.9748 


592.4761 


0.7323 


51 


0.59 


1 


175 


9 


593.9732 


594.4744 


0.7766 


50 


0.77 


1 


180 


10 


604.0442 


604.5450 


0.9996 


49 


0.76 


1 


479 




Fig. 2. Median NACO image of the 51 frames of ^ Vel obtained 
on 2008-04-01. 



With E the number of elapsed orbits. To retrieve the phases listed 
in Table [2l we converted the modified Julian dates of our ob- 
servations into heliocentric Julian dates using the tools by Dan 
Brutor 3- The raw images were processed using the Yoric^EI and 
IRAIosoftware packages in a standard way. A sample image of 
6 Vel is presented in Fig. [2l 

3. Image analysis 

3.1. VISIR 

We obtained photometry of ^Vel A in the PAHl and PAH2 
bands using a circular aperture of 0.3'' in radius, centered on 



^ http://www.physics.sfasu.edu/astro/javascript/hjd.html 
^ http://yorick.sourceforge.net/ 

^ IRAF is distributed by the NOAO, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 
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Fig. 3. Apertures (circles) used for the VISIR photometry of 
^ Vel A (left) and B (right) shown respectively on the PAHl im- 
age #B and the subtracted image #B - #D (see text for details). 



the star position (Fig. [3l left). This small radius is intended to 
avoid the contamination of the photometry of A by component 
B. For ^Vel B, we first subtracted a scaled version of the cali- 
brator star at the position of A from the image before measuring 
the flux through aperture photometry (Fig.[3l right). We thus re- 
moved the contribution from the diflraction pattern (Airy rings) 
of ^ Vel A at the location of B. This process left almost no resid- 
ual thanks to the stability of the PSF. 

We checked that the VISIR images of S Vel and its calibra- 
tors present no detectable difl'use background. This verification 
was achieved by measuring the average flux at distances larger 
than V from the central stars, to avoid the inclusion of the high- 
order Airy rings. Based on this result, we did not apply a sky 
subtraction to our observations, and as we used the same aper- 
tures for the photometry of ^Vel and the flux calibrators, no 
aperture correction was required either. This non-detection of 
the difl'use background was expected for two reasons: firstly due 
to the chopping-nodding observing procedure that removes most 
of the large scale galactic IR background, and secondly because 
the sensitivity of our images is insufficient to detect the spatially 
extended flux from the bow shock in the VISIR bands. The total 
flux from the model proposed by Caspar et al. ( 2008 ) is at most 
3-4 mJy in this wavelength range (Sect. [5]), spread over a surface 
of several arcmin^, far beyond our sensitivity limit. 

The integrated fluxes are afl'ected by difl'erent atmospheric 
absorption due to the difl'erent airmasses compared to the cali- 
brators. To take this into account, we used the empirical formula 
by Schiitz & Sterzik (120051) , that gives the multiplicative correc- 
tion C(A, AM) to apply to the raw flux to remove the atmospheric 
absorption: 



C(A,AM) = 1 -h 



0.220- (^-8.6yum) 



(AM-1) (3) 



We corrected separately the difl'erent observations of S Vel and 
its calibrators. 

In order to absolutely calibrate the measurements, the PAHl 
and PAH2 fluxes from the calibrator HD 67582 were retrieved 
from the spectral template of Cohen et al. (11999D at the central 
wavelength of our filters. For HD 80007 (HIP 45238), we ad- 
justed a model spectrum from Castelli & Kurucz (l2QQ3b to all the 
photometry available in the literature to retrieve its irradiance in 
the Nell filter. We obtained 5.62 Jy (1.02 10"^^ W/m^yum) that 
is within 7% of the value published on the VISIR instrument 
web pag43 (5.26 Jy). The statistical and calibration uncertainties 



http ://w w w. eso. org/sci/f acilities/paranal/instruments/visir/ 
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Table 3. Measured thermal IR irradiances of 6 Vel A and B in 
S.I. units and in Jy. The relative uncertainty cr on the irradiance 
values is listed in the last column. 



Star 


Filter 




10-^^ W/m^/yL/m 


Jy 


cr 


A 


PAHl 


8.59 


36.93 + 2.68 


9.10 + 0.66 


7.3% 


A 


PAH2 


11.25 


12.38 + 1.22 


5.24 + 0.52 


9.9% 


A 


Nell 


12.81 


6.70 + 0.74 


3.68 + 0.41 


11.1% 


B 


PAHl 


8.59 


3.82 + 0.28 


0.94 + 0.07 


7.4% 


B 


PAH2 


11.25 


1.38 + 0.14 


0.58 + 0.06 


10.4% 


B 


Nell 


12.81 


0.56 + 0.19 


0.31+0.11 


34% 



of the measurements were estimated from the dispersion of the 
different available exposures (e.g. #ABEF for the PAHl filter). 
The uncertainty on the Nell measurement, for which only one 
image is available, was taken as the dispersion of the calibrated 
flux for apertures of 150, 300, 600 and 1200 mas. The resulting 
irradiances are presented in Table [3] 

The ratios p = f(B)/f(A) of the fluxes of S Vel B and S Vel A 
in each band can be estimated more accurately than the absolute 
flux of each star due to the removal of the calibration uncertainty. 
We obtain the following values in the three filters: 



P8.6;.m= 10.3 + 0.2%, 
Pll.25/.m= 11.1+0.4%, 
Pl2.81/.m = 8.4±2.3%. 



(4) 
(5) 
(6) 



Fig. 4. Aperture (15 pixels = 0.2''in radius) used for the integra- 
tion of the flux of ^Vel A (left) and B (right), the latter on the 
ring-median subtracted NACO image (see text for details). 
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Fig. 5. Ratio p(2.17) = f(B)/f(A) at 2.17 yum as a function of the 
aperture radius (in pixels) used for the integration of the flux of 
A and B. An aperture of 15 pixels was selected for the analysis 
(square). 



3.2. NACO 

We first obtained classical aperture photometry of ^Vel A us- 
ing circular apertures of 1 to 40 pixels in radius (single pixel 
to 0.5 3 ''radius), for each of our ten epochs. Although the two 
components of ^Vel are well separated in the NACO images, 
the diff'use halo from the residuals of the adaptive optics correc- 
tion of ^ Vel A are not negligible at the position of B. For this 
reason, the photometry of B cannot be measured directly on the 
images. To remove the diff'use wings of the point spread function 
(PSF) of A, we subtracted from each point of the image its ring 
median at the corresponding radius from S Vel A. With this pro- 
cedure, we cleanly subtracted the diff'use background from A, 
mostly without introducing additional noise. Very little residuals 
from star A are present on the subtracted image, thanks to the 
good circularity of the PSF (Fig.|4]). We then measured aperture 
photometry using the same 1-40 pixels apertures as for A, cen- 
tered on B, and computed the flux ratio p = f(B)/f(A) of the 
pair. We do not need to correct for the airmass, as both A and B 
are afl'ected by the same atmospheric absorption. We checked in 
the images that the sky background is negligible due to the short 
exposures and narrow-band filter. 

The choice of the aperture radius is an important factor for 
the accuracy of the flux ratio measurement. In order to obtain the 
best match between the two stars, it is essential to have the same 
integration radius for both stars. Over the 0.6''distance between 
the two stars, the variation of the PSF shape and Strehl ratio 
is negligible, especially for the good seeing conditions of our 
observations. Fig. [5] shows the variation of the derived flux ratio 
as a function of the aperture radius. There is a clear inflexion 
point around 15 pixels in radius (0.2''). This corresponds to the 
radius at which the wings of the PSF of B reach the noise level. 
The increase of p observed for larger apertures corresponds to 
the inclusion of the residuals from A in the photometry. 



To estimate the error bar of the flux ratio, we added quadrati- 
cally the uncertainties due to the choice of aperture radius over a 
10- 15 pixel range (considered as systematic, thus not averaging 
out), and the dispersion of the measurements over the first nine 
epochs (statistical). The flux ratio is assumed to be constant out- 
side of the eclipses, in order to average the nine measurements 
(but see also Sect. 14.41) . For the epoch of the eclipse, we consid- 
ered as uncertainty the standard deviation of the ratios obtained 
at the other epochs. 

From this procedure, we obtain outside of the eclipse: 

P2.i7^m = 9.66 ±0.05% (7) 

corresponding to a magnitude difl'erence of Am = 2.537 ± 0.005. 
During the primary eclipse of ^ Vel A, the flux ratio becomes: 

P2.i7/.m(Eclipse) = 14.50 ± 0.12% (8) 

giving Am(Eclipse) = 2.097 ± 0.009. The eclipse depth at 
2.17 yum is therefore d(Am) = 0.440 ± 0.011 mag. As a remark, 
we did not deconvolve our NACO images, as we did not obtain 
a PSF calibrator for this purpose, and we leave the astrometric 
analysis of these images for a future article. 

4. Photometric analysis 

4.1. Visible magnitudes 

If not specified otherwise, we assume here that the photomet- 
ric measurements of the ^Vel system available in the litera- 
ture always include the A and B components, as their angular 
separation has been smaller than 2" for at least five decades. 
BesseU (.1990 ) gives: 

my (Aab B) = 1 .96 ± 0.02 (9) 
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Fig. 6. Measured value of f(B)/f(A) at 2.17 jim for our 10 mea- 
surement epochs. Note the primary eclipse occuring at the last 
epoch. 



with an associated uncertainty of ±0.02 mag. 

The photometric observations from Hipparcos (ES A 119971) 
give the differential magnitude of Aab (combined) and B in the 
Hp band, that are sufficiently close to V to consider the values 
identical within their error bars: 



Am//p(Aab/B) ^ Amv(Aab/B) = 3.58 ± 0.07 



Considering Bessel's data, this translates immediately into the 
following magnitude for Aab and B: 



my(Aab) = 2.00 ± 0.02 
my(B) = 5.54 + 0.08 



(11) 
(12) 



In addition, the combined visible photometry available in the 
literature (see the compilation by Otero et al. 2000 ) gives us the 
combined Aab+B magnitudes, the depth of the primary (I) and 
secondary (II) eclipses: 



Amy(I, Aab + B) = 0.51+0.01 
Amy(II, Aab + B) = 0.32 + 0.01 



(13) 
(14) 



By removing the flux contribution of B computed precedently, 
we obtain the following eclipse depths for Aab only: 



Amy(I, Aab) = 0.53 + 0.01 
Amy(II, Aab) = 0.33 + 0.01 



(15) 
(16) 



As shown by the curve presented by Otero et al. (12000b , the sec- 
ondary eclipse of d Vel Aab is total, and we can thus derive im- 
mediately the visible magnitudes of Aa and Ab: 

my(Aa) = 2.33 + 0.03 (17) 
my (Ab) = 3.44 + 0.03 (18) 

4.2. Infrared magnitudes 

In addition, our NACO observations (Sect. 13.2b give us the rela- 
tive fluxes of Aab and B in the narrow band filter band at 2. 17 [im 
with a FWHM of 0.023 yum: pi.iitim = 9.66 + 0.05%. As this 
wavelength corresponds to an absorption line of hydrogen that 
is present in the spectra of hot stars, the conversion of the mea- 
sured narrow-band ratios to standard K band ratios requires tak- 
ing the position of the quasi-monochromatic wavelengths within 
the bands and the shape of the observed spectra into account. 

We thus proceeded in two steps to obtain the conversion be- 
tween the narrow-band flux ratio and the broadband flux ratio. 
Starting from the narrow-band fluxes measured on ^ Vel's com- 
ponents, we computed the radius and efl'ective temperature of 



each star assuming (as a first approximation) that the K flux ra- 
tios are identical to the 2.17yum ratios, using the method des- 
cussed in Sect. 14.51 This gave us a first estimation of the spec- 
tral types of the stars. We then used the Pickles (19980 refer- 
ence spectra corresponding to their spectral types (A2IV, A4V 
and F8V) to recover the corresponding broadband flux ratio. 
The K band standard filter profile was taken from Bessell & 
Brett (11988b . We observed that the true ratios in the K band are 
very close to the ratios in the 2.17 yum filter, within 0.01 mag for 
all three stars, which is small compared to the other uncertain- 
ties. We thus neglect this diff'erence in the following analysis. 

The COBE/DIRBE instrument measured Fk = 121.1 ± 
11.4Jy at 2.2 yum (Smith et al. 2004), corresponding to mk = 
1 .77 + 0.10, and the 2MASS catalog (Skrutskie et al. l2QQ6b gives 
an apparent K band magnitude of jtik = 1.72 + 0.26 for Aab+B, 
equivalent to Fk = 127 ± 34 Jy. Averaging these two measure- 
ments, we obtain mk = 1.76 + 0.10 and consequently: 



mK(Asib) = 1.86 + 0.09 



(10) m/r(B) = 4.40 + 0.09 



(19) 
(20) 



4.3. Photometry during the primary eclipse 

The primary eclipse of Aab in the K band observed with NACO 
can give us the diff'erence in surface brightness of Ab and Aa. As 
this eclipse is not total, but annular, we cannot derive unambigu- 
ously the brightness ratio of the two stars as it depends on their 
angular diameters. In this paragraph, we consider that the eclipse 
is perfectly centered, thus neglecting a possible diff'erential eff'ect 
of the limb darkening (LD) of the two stars. As we observed very 
close to the eclipse minimum and in the K band, where the LD 
is small (the LD/UD correction is less than 3%), this appears as 
a reasonable assumption. We also neglect the possible distortion 
of the photospheres of Aa and Ab due to a possible fast rotation. 

During the primary eclipse, the disk of Ab "replaces" part of 
the disk of Aa with a lower eff'ective temperature (hence lower 
surface brightness) "patch". The eclipse depth in the K band 
is 0.44 + 0.01 mag (Sect. 13.2b . This can be compared with the 
depth of the primary eclipse in the visible: Amy = 0.53 + 0.01. 
The "color" of the primary eclipse is therefore Amy - Am^ = 
0.09 + 0.02 mag. This corresponds to the "diff'erential diff'erence" 
in color between the two stars between the visible and IR wave- 
lengths. This leads to a magnitude diff'erence of: 

m^(Ab) - m^(Aa) = 1.02 + 0.04 (21) 

corresponding to a flux ratio /(Ab)//(Aa) = 39.2 + 1.4 % and: 

m^(Aa) = 2.22 + 0.09 (22) 



m^(Ab) = 3.24 + 0.09 



4.4. Emergence of Ab from the secondary eclipse 



(23) 



It is interesting to notice that our 2008-04-23 NACO observa- 
tion has been obtained only about 30 minutes after the end of the 
secondary (total) eclipse during which the 6 Vel Aa component 
covered Ab. The apparent disks of the two stars are therefore 
"touching" each other. Our measurement shows a slight vari- 
ation of the photometric ratio p = f(B)/f(A), as we measure 
9.90 + 0.12%, compared to an average value over the eight other 
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out-of-eclipse epochs of 9.63+0.05%. This means that the mask- 
ing of part of the close environment of d Vel Ab by Aa results in 
a slight loss of total flux of the Aab system. 

Although of marginal statistical significance (2.1 cr), this 
measurement is the maximum value of the out-of-eclipse ratios 
we obtained. It could be due to the presence of a circumstellar 
flux contributor of unknown nature in the immediate vicinity of 
^ Vel Ab, accounting for 0.27 + 0.13% of the total flux of Aab. 
If Ab is a fast rotator, and its rotation axis is perpendicular to the 
orbital plane, one hypothesis could be that it is the signature of 
a circumstellar disk made of gas (Ae episode). 

4.5. Angular diameters, radii and effective temperatures 

From our spatially resolved photometry of the three stars of S Vel 
(Aa, Ab and B), we can deduce their limb-darkened angular di- 
ameters ^LD, photospheric radii R and efl'ective temperatures Tq^ 
using the surface brightness-color (SBC) relations calibrated by 
Kervella et al. (2004) in (V, V - K) and the Hipparcos paral- 
lax. The results are summarized in Table IH This approach al- 
lows us to rely only on empirical SBC relations (calibrated by 
interferometric angular diameter measurements) to compute the 
parameters of the three stars. Another approach to determine 
these parameters would be to adjust directly synthetic spectra 
to the photometry. While the comparison with such model spec- 
tra is necessary to determine the IR excess (Sect. [5]), it has at 
this stage the disadvantage of making the derived parameters de- 
pendent on the choice of one particular numerical model library. 
For this reason, we prefer to use the empirical SBC relations. In 
any case, the good agreement between the model spectra and the 
measured photometry in the visible and near-IR domains shown 
in the bottom parts of Fig. [7] and [8] indicates that both methods 
give consistent results. 

The parameters of 6 Vel are in good agreement with the val- 
ues deduced for the combined system A-fB from spectroscopy 
by Allende Prieto & Lambert ( 1999): T^^ = 8700 K,R = 3.9 Rq. 
The larger radius is naturally explained by the combination of 
the flux from Aa and Ab (and to a lesser extent of B), while the 
efl'ective temperature is a weighted average of the values of the 
three stars. They also find an eff'ective gravity of logg = 3.66 
and a mass of M = 2.51 M©. The reliability of these values is 
uncertain due to the eff'ect of the orbital motion of Aab on the 
spectral line broadening. However, we do not confirm the large 
radii derived by Kellerer et al. (2007) for Aa and Ab (6.0 and 
3.3 Rq) from interferometric observations of the trio. 



5. Discussion 

The IR excess detected around S Vel by IRAS was recently at- 
tributed by Caspar et al. (2008) to interstellar dust heated by 
^ Vel. In this Section, we examine the possibility that one of the 
components of ^Vel A or B presents an IR excess of circum- 
stellar origin. For this purpose, we compare the fluxes extracted 
from our spatially resolved VISIR observations to the expected 
photospheric flux from the Aab and B components. 

To retrieve the synthetic stellar spectra of the three stars from 
the library assembled by Castelli & Kurucz (2003), hereafter 
CK03, we used the parameters listed in Table IH We read the 
tables for log^ values of 4.0 for (5 Vel Aa and Ab, and 4.5 for 
B, and solar metallicities. The sensitivity of the models to small 
variations of these two parameters is minimal. We would like 
to emphasize that we did not fit these spectra to the available 
photometry, but we simply multiplied the CK03 spectra by the 




ii 






• 











Wavelength (mic) 

Fig. 7. Spatially resolved photometry of S Vel Aab and B com- 
pared to Castelli & Kurucz ( 2003 ) models of Aa-hAb and B, 
considering the stellar properties listed in Table |4] (top) and the 
residuals compared the observations (bottom). 

Table 5. Observed fluxes and excesses of S Vel compared to the 
Castelli & Kurucz (2003) models. Ncr is the measured excess 
expressed in number of times the measurement uncertainty. 



Band 


A (fim) 


F (W' W/m'/fim) 


F/F^,d (%) 


Ncr 


Aab+B 
PAHl 
PAH2 
Nell 


8.59 
11.25 
12.81 


4.07 + 0.27 
1.38 + 0.12 
7.26 + 0.08 


99 + 7 
98 + 9 
86+11 


-0.2 
-0.3 
-1.3 


Aab 
PAHl 
PAH2 
Nell 


8.59 
11.25 
12.81 


3.69 + 0.27 
1.24 + 0.12 
0.67 + 0.07 


98 + 7 
97 + 10 
88 + 11 


-0.2 
-0.3 
-1.1 


B 

PAHl 
PAH2 
Nell 


8.59 
11.25 
12.81 


0.382 + 0.028 
0.138 + 0.014 
0.056 + 0.019 


104 + 7 
109 + 10 
74 + 34 


0.5 
0.9 
-0.7 



squared limb-darkened angular diameter of each star (Table |4l). 
From these spectra, we derived the expected thermal IR flux 
from the stellar photospheres of Aa-i-Ab and B in the PAHl, 
PAH2 and Nell filters corresponding to our VISIR observations. 

We can also add the CK03 spectra of the three stars in or- 
der to compare the predicted flux of the trio to the existing 
non spatially resolved measurements of A-hB from the litera- 
ture. The resulting spectrum is shown in Fig. [SI The photomet- 
ric measurements are taken from Thompson et al. (|1978t UV 
fluxes from TDl), Morel et al. ([T978t UBVRI), Skrutskie et 
al. (2006; 2MASS JHK), Smith et al. (2004; COBE/DIRBE, 
four wavelengths from 1.25 to 4.5 yum), VISIR (8.6, 11.25 and 
12.81 yum, this work), Jayawardhana et al. (120011 N and 18.2 yum 
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Table 4. Properties of d Vel Aa, Ab and B derived from their V and K band magnitudes using the surface brightness-color relations 
from Kervella et al. (12004b . The bolometric corrections for the V band were taken from Bessell et al. (11998b . The approximate 
masses of the three stars are discussed in Sect. [51 



Star 


^LD (mas) 


Radius (Rq) 


reff(K) 


M 


V 


BCy 


Mboi 




Spect. 


M/Mq 


Aa 
Ab 
B 


1.21 + 0.03 
0.77 + 0.02 
0.53 + 0.02 


3.17 + 0.08 
2.02 + 0.05 
1.39 + 0.06 


9000 + 400 
8 600 + 350 
6350 + 350 


0.39 + 
1.50 + 
3.60 + 


0.04 
0.04 
0.04 


-0.07 
0.00 
0.00 


0.32 + 0.04 
1.50 + 0.04 
3.60 + 0.04 


59.2 + 2.0 
20.0 + 0.7 
3.1+0.2 


A2IV 
A4V 
F8V 


2.5 + 0.1 
2.0 + 0.1 
1.3 + 0.1 





3,95 3,90 3,85 

Log Teff (K) 



3,75 3,70 



200% 
180% 
o 160% 
£ 140% 

LL 

o 120% 
ir 100% 









T / T 






r : T / t 












It in 




y 






a. ■*> — cr 


[ 















Wavelength (mic) 

Fig. 8. Total model flux of S Vel Aab-hB (solid curve), together 
with the available photometry from the literature (squares) and 
the VISIR photometry (filled circles). The dashed curve is the 
flu x contr ibution of the IR excess model proposed by Caspar et 
al. C2008'). The bottom plot shows the relative photometry with 
respect to the combined photospheric model of the three stars. 



bands), Caspar et al. (l2QQ8l Spitzer/MIFS 24 and 70yum), and 
IPAC ( 1986; IRAS 12, 25 and 60 /im). The additional flux from 
the bow shock model proposed by Caspar et al. ( 2008) is shown 
on both panels of Fig. [8] as a dashed curve. Its contribution in 
the range of wavelengths sampled by our VISIR observations 
appears negligible (^0.1% of the photospheric flux). 

As shown in Fig. [71 Fig. [8] and Table [5j the comparison 
with our VISIR spatially resolved photometry of S Vel A and B 
does not show a significant circumstellar IR excess in the VISIR 
bands compared to the expected photospheric flux of each of the 
two components, at a few percent level. The excess that is visible 
long wards of lOyum in Fig. [8] comes from the flux contribution 
of the ISM bow shock at large angular distances that has been 
observed by Caspar et al. (^2008j . We note that the IRAS and 
Jayawardhana et al. (120011) photometry shows an excess com- 
pared to their model between 12 and 25 jum. This may be a con- 
sequence of the fact that the flux contribution from the ISM is (at 
least partly) included in the large aperture measurements (IRAS 
in particular), but not in our narrow aperture photometry. 



Fig. 9. Positions of S Vel Aa, Ab and B in the HR diagram, with 
the isochrones and isomass curves from Cirardi et al. (12000b 
overplotted. 

We notice a relatively large ultraviolet excess on ^Vel in 
Fig. [SI A fast rotating star usually presents a UV excess due to its 
overheated polar caps (von Zeipel efl'ect, see e.g. Aufdenberg et 
al. 120061) . A discussion on the rotational velocity of the compo- 
nents of S Vel A can be found in Argyle, Alzner & Horch (12002b . 
As a remark, Royer, Zorec & Comez (12007b found a rotational 
velocity of vsin/ = 150km/s, but this was computed for the 
combined spectrum of the three stars of S Vel so a bias may be 
present. 

The physical parameters derived for the three stars (Table |4l) 
indicate that ^Vel Aa already evolved away from the main 
sequence. This is confirmed by the evolutionary tracks from 
Cirardi et al. (2000), shown in Fig. [9l that give approximate 
masses of 2.5, 2.0 and 1.3 M© respectively for Aa, Ab and B, and 
an age of approximately 400 - 500 Myrs for the system. This age 
compares well with the age of 330-390 Myr proposed by Rieke 
et al. (12005b . while the masses are in relatively good agreement 
with Argyle, Alzner & Horch ('2002'), who determined a total dy- 
namical mass of the 6 Vel AB system of 5.7^| j M©. We postpone 
a more detailed discussion of the ages and evolutionary status of 
the S Vel stars to a forthcoming paper. 

6. Conclusion 

Our VISIR photometry of S Vel A and B does not show the pres- 
ence of an excess of circumstellar origin in the thermal IR do- 
main (8-13 Jim) at a level of a few percents. This result indicates 
that these stars do not host a warm debris disk with a typical 
temperature around 200-300 K. This supports the conclusions of 
Caspar et al. (2008) who attribute the IR excess detected with 
Spitzer to the surrounding interstellar medium material. The ab- 
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sence of a warm circumstellar disk in the inner d Vel system may 
be a consequence of the gravitational interactions between the 
three stars. They could have caused the dispersion of the residual 
circumstellar material from which they stars formed. However, 
the possibility still exists that a cold debris disk is present, as in 
the case of Vega (Su et al. 12005b and Fomalhaut (Stapelfeldt et 
al. 2004), but at large distances from the stars. In order to test 
this scenario, observations of ^ Vel in the far IR or submillimet- 
ric domain similar for instance to those of Fomalhaut by Holland 
et al. (2003 ) would be necessary. 

From the flux ratios in the V and K bands, we could derive 
the physical properties of the three stars Aa, Ab and B. ^Vel 
appears as a moderately evolved system, with the Aa compo- 
nent currently leaving the main sequence. As suggested by the 
V sin / value of Royer et al. (I2007L the ultraviolet excess of 6 Vel 
and our putative observation of near-IR circumstellar emission 
close to Ab may indicate that 6 Vel Ab (and/or Aa) could be a 
fast rotator, turning this eclipsing system into a promising object 
to map the polar brightening associated to the von Zeipel efl'ect 
(von Zeipel|1924i). 

Acknowledgements. We are grateful to the ESO observing team at Paranal for 
the perfect execution of our NACO observations of (5 Vel, in particular for the 
time-critical eclipse on May 18th, 2008. We also warmly thank Sebastian Otero 
(Asociacion Cielo Sur, Argentina), Brian Fraser (Sunninghill Observatory, South 
Africa) and Jaime Garcia (Instituto Copernico, Argentina) for providing the 
visible photometry of (5 Vel. This work received the support of PHASE, the 
high angular resolution partnership between ONERA, Observatoire de Paris, 
CNRS and University Denis Diderot Paris 7. This research took advantage of the 
SIMBAD and VIZIER databases at the CDS, Strasbourg (France), and NASA's 
Astrophysics Data System Bibliographic Services. 



Skrutskie, R. M., Cutri, R., Stiening, M. D., et al. 2006, AJ, 131, 1163 
Smith B. J., Price S. D., Baker R. I. 2004, ApJ Suppl. Ser, 154, 673 
Stapelfeldt, K. R., Holmes, E. K., Chen, C, et al. 2004, ApJS, 154, 458 
Su, K. Y. L., Rieke, G. H., Misselt, K. A., et al. 2005, ApJ, 628, 487 
Thompson G. I., Nandy K., Jamar C, et al. 1978, "Catalogue of stellar ultraviolet 

fluxes (TDl)", The Science Research Council, U.K. (1978) 
Tokovinin, A., Sarazin, M. & Smette, A. 2007, MNRAS, 378, 701 
von Zeipel, H. 1924, MNRAS 84, 665 



References 

Allende Prieto, & Lambert D. L. 1999, A&A, 352, 555 

Argyle, R. W., Alzner, A., & Horch, E. R 2002, A&A, 384, 171 

Aufdenberg, J. P., Merand, A., Coude du Foresto, V, et al. 2006a, ApJ, 645, 664, 

Erratum 2006, ApJ, 651, 617 
Bessell, M. S. & Brett, J. M. 1988, PASP, 100, 1134 
Bessell M. S., Castelli F., & Plez B. 1998, A&A, 333, 231 
Bessell, M. S. 1990, A&A Suppl. Ser, 83, 357 

Castelli, F. & Kurucz, R. L. 2003, lAU Symposium 210, Modelling of Stellar 

Atmospheres, Uppsala, Sweden, eds. N.E. Piskunov, W.W. Weiss, and D.F. 

Gray, 2003, ASP-S210, arXiv:astro-ph/0405087 
Cohen, M., Walker, R. G., Carter, B., et al. 1999, AJ, 117, 1864 
Doucet, C, Lagage, R O. & Pantin, E. 2007, Proc. of the VIRA conference, 

Paris, 20-22 march 2006, in press, astro-ph/06 10322 
Domiciano de Souza, A., Kervella, P., Bendjoya, Ph, & Niccolini, G. 2008, 

A&A, 480, L29 
Doucet, C, Habart, E., Pantin, E. et al. 2007, A&A, 470, 625 
ESA 1997, The Hipparcos and Tycho Catalogues, ESA SP-1200 
Caspar, A., Su, K. Y. L., Rieke, G. H., et al. 2008, ApJ, 672, 974 
Girardi, L., Bressan, A., Bertelli, G., Chiosi. C. 2000, A&A Suppl. Ser, 141, 371 
Holland, W. S., Greaves, J. S., Dent, W. R. F., et al. 2003, ApJ, 582, 1141 
IPAC 1986, IRAS Catalog of Point Sources, Version 2.0, Joint IRAS Science 

W.G. 

Jayawardhana, R., Scott Fisher, R., Telesco, C. M., et al. 2001, AJ, 122, 2047 
Kellerer, A., Petr-Gotzens, M., Kervella, P., & Coude du Foresto, P. 2007, A&A, 
469, 633 

Kervella, R, Thevenin R, Di Folco, E., & Segransan, D. 2004, A&A, 426, 297 
KerveUa, R, & Domiciano de Souza, A. 2007, A&A, 474, L49 
Lagage, RO. et al. 2004, The ESO Messenger 117, 12 

Lenzen, R., Hofmann, R., Bizenberger, R, & Tusche, A. 1998, SPIE 3354, 606 
Masciadri, E., Brandner, W., Bouy, H. et al. 2003, A&A, 411, 157 
Morel M., & Magnenat R 1978, A&A Suppl, 34, 477 
Otero, S. A., et al. 2000, IBVS, 4999 
Pickles, A. J. 1998, PASP, 110, 863 

Rieke G. H., Su K. Y. L., Stansberry J. A., et al. 2005, ApJ, 620, 1010 

Royer, F., Zorec J., & Gomez A. E. 2007, A&A 463, 671 

Rousset, G., Lacombe, F., Puget, F., et al. 2003, Proc. SPIE 4839, 140 

Schiitz, O., & Sterzik, M. 2005, Proc. ESO Workshop "High Resolution Infrared 

Spectroscopy in Astronomy", H. U. Kaufl, R. Siebenmorgen, and A. 

Moorwood ed., 104-108 



